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Systems evolving under the influence of competing two-body and three-body interactions, are of 
particular interest in exploring the stability of eqnilibrium states of strongly interacting many-body 
system. We present a solvable model based on qubit networks, which allows us to investigate the 
intricate influence of these couplings on the possible asymptotic equilibrium states. We study the 
asymptotic evolution of finite qubit networks under two and three-qubit interactions. As represen¬ 
tatives of three-qubit interactions we choose controlled unitary interactions (cu-interactions) with 
one and two control qubits. It is shown that networks with purely three-qubit interactions exhibit 
different asymptotic dynamics depending on whether we deal with interactions controlled by one 
or two qubits. However, when we allow three-qubit interactions next to two-qubit interactions, the 
asymptotics is dictated by two-qubit interactions only. Finally, we prove that the simultaneous 
presence of two types of three-qubit interactions results in the asymptotic dynamics characteristic 
for two-qubit cu-interactions. 

PACS numbers: 03.65.Aa, 03.67.-a 


I. INTRODUCTION 

The evolution of any physical system is governed by the 
interactions between the constituting parts and the ex¬ 
ternal conditions. Especially in statistical physics the un¬ 
derstanding of the form of the internal interactions plays 
a crucial role in the description of the system. Usually 
we limit ourselves to two-body interactions and neglect 
higher order (many-body) interactions. The ignored in¬ 
teractions have usually properties which differ from two- 
body interactions and induce different behavior of the 
system. In particular, states which are stable under two- 
body interaction may become unstable under the influ¬ 
ence of additional three-body interaction as observable in 
undercooled gases, for example. Hence there is certainly 
interest in situations in which competing types of inter¬ 
actions can be solved analytically and allow to specify 
the role of the additional contributions. We attack this 
problem by using the framework of quantum networks 
and study the dynamics of arbitrary sized network with 
competing two-body and three-body interactions. The 
elementary interactions are given in terms of quantum 
gates. 

Qubit systems play an important role in quantum the¬ 
ory. Beside their importance as nontrivial representa¬ 
tions of many physical systems they are of significant use 
in quantum information theory m- Coupled qubits 
equipped with a special type of interaction can be used 
in quantum computing to process and manipulate infor¬ 
mation. The mutual couplings do not necessarily involve 
only particular qubits, they can be also a result of an 
interaction with an additional auxiliary system. Qubit 
systems of this kind are generalizations of classical net¬ 
works HIS]. They are therefore called qubit networks and 
they can be represented by graph structures. Nodes are 
qubits (equipped with two base states |0) and |1)), edges 
(directed or undirected) represent interactions between 


nodes. Qubit networks are a part of a broader family of 
quantum networks. Because of their importance, quan¬ 
tum networks received a lot of attention in recent years 
[SHE]- They offer a wide variety of utilizations, e.g. es¬ 
tablishing entanglement for various applications HE], 
modelling of quantum decoherence [laiis] or reaching 
quantum consensus [IlHISI- 

We often find ourselves in situation, in which we are 
unable to control the type of interaction or the strength 
of the coupling between qubits. As a result, classical ran¬ 
domness must be introduced into such a qubit network. 
A convenient way how to describe such a dynamics are 
Markov chains, i.e. sequences of random (quantum) oper¬ 
ations. In a significant number of cases the situation can 
be mathematically described by a superoperator which 
represents a random unitary operation (RUO) [17j . This 
superoperator is equipped with a probability distribution 
over all possible qubit couplings that can be realized in a 
single application of the given RUO. The resulting time 
evolution is then given by an iteration process of repeated 
application of the given RUO on the density matrix de¬ 
scribing the state of the system in the previous step. 

Recently, quite a general class of two-qubit interactions 
has been considered. The interaction between qubits has 
been chosen as a member of a one-parameter family of 
two-qubit cu-interactions (i.e. controlled unitary interac¬ 
tions) OIIH]. The evolution of qubit networks with such 
interaction was calculated in the asymptotic regime, i.e. 
after a sufficiently large number of interactions. The de¬ 
scription had taken into account the topology of the given 
qubit network, i.e. realizable couplings between qubits. 
Such qubit networks can be used to model a number of 
physically relevant and interesting situations, e.g. they 
can be used as a simple model of particle collisions in a 
rarehed gas or as a class of one-parameter decoherence 
models [12]. Since the asymptotic state is independent 
of the non-zero values of the probability distribution over 
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the qubit couplings it is of interest to find out, how 
does the asymptotic state change, if we enrich the dynam¬ 
ics of the network by adding higher-order interactions, 
e.g. three-qubit interactions as the simplest generaliza¬ 
tion. One is also interested in properties of interactions 
of higher orders on their own, i.e. without the presence of 
two-qubit interactions, as the results obtained previously 
cannot be extended to more general cases in a straight¬ 
forward way. Higher order interactions can be considered 
a step towards more complicated interaction structures 
which are not decomposible into simple (random) two- 
qubit interaction. 

In the following, we investigate the effect of three-qubit 
cu-interactions on the qubit network equipped with two- 
qubit cu-interactions. The considered three-qubit cu- 
interactions have form of controlled unitary three-qubit 
interactions with one or two control qubits. Both cases 
are studied in detail for a class of interaction topologies 
represented by so-called base graphs. First, we generalize 
the concept of an interaction graph [T^ (used to spec¬ 
ify realizable two-qubit couplings) to the case of three- 
qubit couplings. Next, we derive the form of asymptotic 
states of three-qubit cu-interactions for certain family of 
interaction topologies. This enables us to study proper¬ 
ties of qubit networks equipped with purely three-qubit 
cu-interactions as well as properties of qubit networks 
equipped with both two and three-qubit cu-interactions. 

The paper is structured as follows. Section pro¬ 
vides the theoretical background of basic concepts of 
graphs and hypergraphs, which are used through the 
rest of the paper. Section im] summarizes basic prop¬ 
erties of RUOs which are required for finding the asymp¬ 
totic state (for more details see [HI)- Results relevant to 
two-qubit cu-interactions are given in section EYl Rigor¬ 
ous definitions of three-qubit cu-interactions and corre¬ 
sponding graph structures used through the paper are 
given in section |V] as well as the derivation of corre¬ 
sponding asymptotic states. The relation between two 


and three-qubit cu-interactions is studied in section VI 


Section m presents additional properties of three-qubit 
cu-interactions. Summary and outlook of the paper are 
left to section [Vnl Main steps of proofs of theorems pre¬ 
sented in section are given in appendix. 


II. BASIC DEFINITIONS IN GRAPH AND 
HYPERGRAPH THEORY 


Let us first introduce graph structures which are con¬ 
venient for the description of multi-qubit cu-interactions. 
A natural and often used way to inscribe a topology of 
two-qubit cu-interactions in qubit networks constitute di¬ 
rected graphs. A directed graph is defined as an ordered 
pair G = (V, E), where V = {ui,... utv} is the set of ver¬ 
tices and E = {ei... Cm} is the set of edges. Each edge 
can be expressed as a pair of vertices, which are called the 
tail and the head of the given edge, hence 
A sequence of vertices (u ^^,... ,Vi^) such that (vt ^, ) 


forms an edge in G is called a directed path connecting 
vertices Vi^ and Vi^ . These vertices are called the begin¬ 
ning and the end of the path respectively. If every two 
vertices of the directed graph G are connected by a path 
in both directions G is said to be strongly connected. 

A hypergraph m is a direct generalization of a graph 
in which we do not restrict relations between vertices to 
be binary ones. Instead the directed hypergraph is de¬ 
fined as an ordered pair E[ = {V, E) where E is an arbi¬ 
trary subset of the set constituted by all pairs of disjoint 
subsets X d V. Thus every e G E can be written as 
e = (A, Y) with A, Y C V, A, Y 0 and A n Y = 0. If 
u G X, then u is called the tail of hyperedge e, if u G Y, 
then V is called the head of hyperedge e. A directed 
path in a hypergraph is defined as a sequence of vertices 
(uii, ..., Vij^) such that for every j there is (A, Y) G E 
such that Vj G X and vj+i G Y. A directed hypergraph 
H is called strongly connected if every pair of vertices can 
be connected by a directed path in both directions. 

There exist special types of hyperedges which we find 
particularly useful for describing the topology of multi¬ 
qubit cu-interactions. Let e = (A, Y) be a hyperedge. If 
|A| = 1, e is called a F-arc, if |Y| = I, e is called a B-arc. 
A hypergraph E[ in which all hyperedges are F-arcs is 
called a F-graph, a hypergraph El in which all hyperedges 
are B-arcs is called a B-graph. Ordinary graphs belong 
to both these sets. 

Natural partial ordering on hypergraphs is defined in 
the following way. Let G = {V,E) and H = {U,F) be 
two hypergraphs. We say that G is a subgraph of iL, 
i.e. G C H, a G and E[ are defined on the same set of 
vertices and e G E implies e G F. 


HI. ATTRACTOR METHOD FOR DYNAMICS 
GENERATED BY RUO 

In this paper we study quantum systems undergoing 
discrete dynamics whose one step of evolution is de¬ 
scribed by a random unitary operation (RUO). General 
properties of a random unitary evolution with special em¬ 
phasis on its asymptotic regime was extensively studied 
in m- We briefly summarize main findings important 
for our considerations to follow. 

Assume a quantum system associated with a finite 
dimensional Hilbert space JY and let us denote 
the Hilbert space of all operators acting on the Hilbert 
space Jif (equipped with Hilbert-Schmidt scalar prod¬ 
uct [T]). One step of evolution is given by the RUO 
$ : B{M’) -G H(JY) whose action on the system ini¬ 
tially prepared in a general mixed state p can be written 
in the form 

n 

<^>{p) = Y,p.U,pUl ( 1 ) 

i-1 

with unitary operators Ui and probability distribution 
{Pi}7=i- R-UO <I> is a quantum operation with Kraus op¬ 
erators defined by Ki = . It belongs to the class of 
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trace-preserving unital quantum operations leaving the 
maximally mixed state invariant. From a physical point 
of view, RUO takes into account our lack of information 
which particular unitary evolution the system undergoes 
and incorporates in an incoherent manner all exclusive 
unitary paths of evolution represented by different uni¬ 
tary operators Ui properly weighted with probabilities 
associated with all these paths. 

The evolution of the system results from repeated ap¬ 
plications of RUO $. Starting from the initial state 
p(0), the n-th step of the iterated dynamics reads p{n) = 
^{p{n — 1)). In general, the operator $ is neither a her- 
mitian nor normal linear map and consequently a diag- 
onalization of RUO $ in some orthonormal basis is not 
guaranteed. Fortunately, the latter does not apply to 
the asymptotic part of evolution and one can exploit the 
fact that the asymptotic regime of the evolution takes 
place in the so-called attractor space Atr($) C B{J^) 
constructed as 

Atr($) = 0 Ker($-A/). (2) 

AGcr|i[ 

The attractor spectrum f7|i| denotes the set of all eigen- 
values A of with |A| = 1. For a given A S (T|i| the 
corresponding kernel Ker($ — XI) is formed by the solu¬ 
tion of attractor equations 

UaxjU} = XXx,j Vi. (3) 

Any state in the asymptotic regime takes the form 


p{n > 1) = 


Dx 

E 


A”Tr 


p(o)4,. 


(4) 


Ae<7|i| ,2=1 


where Dx = dim[Ker(<i)—A/)] and {Xx^i\i S {1,... ,Dx}} 
forms an orthonormal basis of the subspace Ker($ — XI) 
with respect to the Hilbert-Schmidt scalar product. Ap¬ 
parently, attractors solving equations (|^ are not affected 
by the particular form of the probability distribution {pi} 
provided that ^ 0. As a direct consequence, the 
asymptotics of given RUO is independent on particular 
values Pi as long as they are different from zero. 

At this point we have to stress that attractors are not, 
in general, density operators, i.e. they do not repre¬ 
sent states. This can cause difficulties in the analysis of 
asymptotic dynamics like fixed points, decoherence-free 
subspaces. To overcome this obstacle one may employ 
the so-called pure-state method [50] . Without going into 
details we recapitulate its main points. Let us denote 
{\4'a,ja)} the orthonormal basis of common eigenstates 
of unitaries {Ui}, i.e. 

(h) 


satisfies equations 

U,XU] = XX Wi,j (7) 

and consequently belongs to the subspace of attractors 
corresponding to the eigenvalue A of RUO $. On the 
other hand, any operator satisfying can be decom¬ 
posed into common eigenvectors as (|^. Attractors which 
can be constructed from common eigenvectors are called 
p-attractors. As they satisfy more restricting set of equa¬ 
tions 0 they do not constitute the whole attractor space. 
In particular, the identity operator is an attractor but 
not a p-attractor (except the trivial case of an unitary 
evolution). The space of attractors always contains, as 
the minimal subspace, the span of p-attractors and the 
identity operator. Surprisingly, in some nontrivial cases 
this minimal subspace forms the whole attractor set and 
the asymptotic dynamics can be analyzed easier. Indeed, 
assume the orthogonal projection V onto the subspace of 
common eigenstates of unitaries {Ui}. Let V be its or¬ 
thogonal complement projection satisfying V + V = I. 
Both projections are fixed points of quantum operation 
Q, they both reduce the random unitary operation Q 
and the asymptotic dynamics of the initial state p(0) can 
be written as 


pin^l) = UrVp{0)V (U^y 


Tr 

pio)r 

Tr 

V 



( 8 ) 


for any pair of indices i,j. In this case the asymptotic 
evolution can be understood as an incoherent mixture 
of the unitary dynamics inside the subspace of common 
eigenstates and the maximally mixed state living on the 
orthogonal complement of this subspace. We often en¬ 
counter the special case cr|i| = {!}, which further simpli¬ 
fies the asymptotic evolution given by (|^ to the form 


Poo 


Vp{0)V -f 


Tr 

p{0)v 

Tr 

V 



(9) 


In such a situation any initial quantum state p(0) evolves 
towards the stationary state (|^. 

Let us hnally point out a crucial property of general 
attractor equations ([^ which we extensively use through¬ 
out the paper. We formulate this property in connection 
with partial ordering on hypergraphs. In order to pro¬ 
ceed consider a hypergraph H = (U, E) and assume that 
for any hyperedge e G E there is an unitary operation 
Ue and nonzero probability Pe such that = 1- 

RUO associated with the given hypergraph is defined as 


where index takes into account the degeneracy of 
common eigenvalue a. Any matrix from the span of 
{\4>a,jc,){4‘0,jfi\} with hxed product a/3 = A, i.e. 

E (6) 

aP=X,jc,Pf3 


e&E 

If two hypergraphs G and H satisfy G C H then attrac¬ 
tor spaces of associated RUOs follow relation Atr($//) C 
Atr($G)- This is an immediate consequence of the fact 
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that the hypergraph H implies more constraints than 
the hypergraph G. These relations imply the existence 
of the set of so-called base attractors. Indeed, assume 
a subset of hypergraphs. Any hypergraph from this set 
determines corresponding RUO (10). The subspace of 
base attractors is identified as the maximal set of attrac¬ 
tors included in attractor spaces of all these associated 
RUOs. This subspace always contains the identity op¬ 
erator. Similarly, we define the set of base hypergraphs 
whose associated attractor spaces are formed from base 
attractors solely. In other words, the attractor subspace 
of all RUOs associated with base graphs is constituted 
from base attractors. 


It is important to stress that the definition of base at¬ 
tractors as well as the definition of base graphs is closely 
linked to the chosen set of hypergraphs. In particular, 
in this paper we always consider the set of graphs or a 
certain subset of F-graphs. 


IV. ASYMPTOTIC DYNAMICS OF 
TWO-QUBIT CU-INTERACTIONS 


where 


= cos{(j))Zj + sin{(t>)Xj (12) 

with the identity operator Ij and Pauli operators Xj and 
Zj acting on the j-th qubit [T] . The parameter (j) is cho¬ 
sen from the interval (0,7r). As these couplings have the 
form of asymmetric control transformations, the interac¬ 
tion graph becomes a directed graph in which the trans¬ 
formation is represented by the directed edge e with 
the vertex i as the tail of e and the vertex j as the head of 
e. Once the interaction graph is given the corresponding 
RUO describing one step of evolution is defined as 




(13) 


eGE 


The asymptotic behavior of dynamics generated by (13) 


can be summarized as follows (for details see [IH])- It 
was shown that base graphs are formed by all strongly 
connected graphs. The space of common eigenvectors for 
strongly connected networks is two-dimensional and its 
basis reads 


Before discussing three-qubit cu-interactions we 
shortly review known results on asymptotic dynamics 
generated by two-qubit cu-interactions which have been 
recently studied nuimiiH]. For this purpose let us con¬ 
sider a qubit network consisting of N qubits in which 
pairs of qubits are randomly coupled to each other via 
certain transformations according to the prescribed prob¬ 
ability distribution. Information about qubit couplings 
can be encoded into a directed graph G = {V, E) which 
is called the interaction graph. Each vertex of this graph 
represents a single qubit, e.g. the vertex i represents the 
qubit i. Edges of the interaction graph then represent 
couplings between qubits, e.g. the edge e = (i,j) rep¬ 
resents the coupling between qubits i and j. This edge 
is present in the interaction graph if only if the corre¬ 
sponding probability of realization of given coupling is 
nonzero. The interaction graph defined above gives us 
sufficient amount of information about the given RUO 
as the precise form of the probability distribution {pi} 
is irrelevant for the asymptotic dynamics generated by 
RUO. 

In view of equations (§, the topology of the inter¬ 
action graph determines the attractor space associated 
with a given two-qubit cu-interaction and consequently 
the asymptotic evolution of any initial state Q- Assum¬ 
ing all interaction graphs the main task is to identify their 
corresponding subsp ace of base attractors and the set of 
base graphs (section |nT| . Previous work was focused on 
a one-parameter family of two-qubit cu-interactions in¬ 
cluding the controlled-NOT transformation as its special 
case. In the computational basis the coupling between 
qubits i and j reads 

= | 0 z ) ( 0.1 0 + |I.) ( 1.1 0 uf, (11) 


\0n) = \0f^, 

= (^cos| |0)-f sin| |l)j 

Both these common eigenvectors correspond to the eigen¬ 
value a = 1. According to section [Inl p-attractors consti¬ 
tute 4-dimensional space and in order to complete the at¬ 
tractor space one has to add the identity operator. Only 
for V = 2 qubits there is also one extra non-p-attractor 
corresponding to eigenvalue A = — 1, namely the operator 
A_i = cos(<()/2)(|01)(ll| - |10)(11|) -f sin(<^/2)|01)(10| + 
-I- h.c. Apart of the special case N = 2, the asymptotic 
dynamics is stationary and takes the form (§. 

The main goal of this paper is to analyze how an in¬ 
terplay among two-qubit and three-qubit cu-interactions 
reshapes the explicit form of base attractors and base 
graphs and determines the resulting asymptotic dynam¬ 
ics. To gain insight into this issue one has to first re¬ 
veal the structure of base attractors and base graphs for 
purely three-qubit cu-interactions. 

V. ASYMPTOTIC DYNAMICS OF 
THREE-QUBIT CU-INTERACTIONS 

The three-qubit cu-interaction does not reduce simply 
to the two-qubit cu-interaction and hence this situation 
requires separate attention. For this reason we exam¬ 
ine three-qubit cu-interactions of two different forms by 
changing the number of control qubits. The main pur¬ 
pose is twofold. First, we want to compare properties of 
asymptotic states of two-qubit cu-interactions with prop¬ 
erties of asymptotic states of three-qubit cu-interactions. 
Second, we intend to reveal the relationship between 
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FIG. 1: An example of an interaction Fi-graph G correspond- 
ing to the RUO ' acting on 4 qubits. It involves three 
unitary operations corresponding to three hyperedges (1, 23), 
(3,14), and (2,13). This Fi-graph is not strongly connected 
and thus does not belong to the set of base graphs. 


base attractor spaces of two-qubit and three-qubit cu- 
interaction networks. In particular, it is an important 
question, whether the base attractor space of two-qubit 
cu-interactions is a subspace of the base attractor space 
of three-qubit cu-interactions. In the case of a positive 
answer, the set of attractor equations ([^ implies that 
by enabling three-qubit cu-interactions next to two-qubit 
cu-interactions, one does not change the asymptotic state 
at all as long as the interaction graph corresponding to 
the two-qubit cu-interaction remains strongly connected. 


A. Controlled unitary three-qubit transformation 
with one control qubit 

Let us consider a qubit network in which qubits are 
coupled to each other by three-qubit transformations 
which in the computational basis take the form 

= |0*) (0.1 ® l3k + |1.) (1.1 ® ® (14) 

with defined by (j^. The set of all assumed three- 
qubit controlled interactions ( |14[ ) is prescribed by an in¬ 
teraction Fi-graph G = {V, section [nT| ). The set of 
vertices V = {1,..., N} corresponds to the set of qubits 
and F-arcs {X,Y) G E with |F| = 2 represent possible 
three-qubit cu-interactions. In particular, if Pi^jk > 0 
then the vertices i,j and k are connected in G by a F-arc 
e G E whose head are qubits j and k and whose tail is 
the qubit i. Hyperedges of any Fi-graph thus consist of 
two heads and one tail. An example of a network rep¬ 
resenting a topology of three-qubit cu-interactions with 
one control qubit is depicted in Fig{^ 

Once an interaction Fi-graph G = {V,E) is given the 
associated RUO describing one step of evolution reads 

eeE 


The eigenvalues of the transformation (14) are a = ±1 
and thus the attractor spectrum (j\i\ of RUO (15) is a 
subset of {—1,1}. The following theorem fully describes 
base Fi-graphs and their corresponding base attractors 
of RUOs (15) associated with interaction Fi-graphs. The 
sketch of the proof is given in appendix. 


Theorem V.l Base Fi-graphs of the considered family 
of interactions are all strongly connected interaction Fi- 
graphs G = (U, E). The space of common eigenvectors is 
three-dimensional with the basis 


|0w) = 10)®"^, 

= (^cos^ |0)+sin^ |1)^ 

Wn) = (^sin^ |0)-cos^ |1)^ 

All these common eigenvectors correspond to the eigen¬ 
value a = 1. Base p-attractors form a 9-dimensional 
subspace of the base attractor space corresponding to the 
eigenvalue A = 1. In order to complete the attractor 
space, several cases need to be distinguished. For (j) ^ 
the base attractor space consists of the base p-attractors 
and the identity operator. If (j) = ^ and N > 3 there are, 
besides p-attractors, two additional linearly independent 
attractors corresponding to eigenvalue A = 1 

In,e= E (l + (-l)"^^^)l^)(^|, 

zei(N) 

In,o= E (i-(-i)^^^^)|^)(^I, 

zei(N) 


with I{N) denoting the set of all possible binary N-tuples 
and t(z) is the sum of the bit values of the N-qubit string 
z. Only in the special case d = f o.iT'd N = 3 qubits 
there is an additional attractor X_i corresponding to the 
eigenvalue A = — 1 given by 

X_i = C^(|101) (0111 - 1110) (011| -k |110) (101|) + h.c. 
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Using theorem |V.1[ one can write down asymptotic 
evolution generated by RUO (15). In the case N > 
3 A d ^ the initial quantum state p(0) approaches 
the stationary state ^ with 7^04) being projector on 
the space of common eigenvectors spanned by vectors 
|0jv), and 


B. Controlled unitary three-qubit transformation 
with two control qubits 

In the second case we consider a qubit network in which 
qubits are coupled to each other by three-qubit transfor¬ 
mations. The interaction (in the computational basis) 
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takes the form 


ultl = ihj - |ia,) (ia,|) ®h + |i*i,) ( 1 * 1,1 0 ur. 

(16) 

with dehned by (12). The considered RUOs have 
the form 


{ 4 >) 




eGE 


(17) 


The family of unitary transformations (16) includes also 
the Toffoli gate. Similar to the previous case, E is the 
subset of the set of partially ordered triples {i,j-,k) = 
{j,i] k). The eigenvalues of the transformation 


ij,k 


are 


again given by a = ±1 and the attractor spectrum a\i\ 
thus fulHls the relation crpi C { — 1,1}. 

B-graphs seem to be a natural choice to describe the 
topology of three-qubit cu-interactions (16). However in 


order to fully describe the topology of base graphs in the 
case of controlled unitary interactions with two control 
qubits, we must significantly change the graph represen¬ 
tations of given interactions. The required graph rep¬ 
resentation should provide crucial information whether 
any pair of control qubits can ’’see” all other qubits via 
interaction with them. This interaction does not have to 
be the direct one, it can be mediated by other pairs of 
qubits. This ’’visibility” is mathematically represented 
by paths in given graph representation. For this reason 
the representation constructed analogically as the pre¬ 
viously used representations (interaction graphs and in¬ 
teraction Fi-graphs) is not sufficient, as paths in such 
hypergraph ignore information about one of vertices cor¬ 
responding to one of control qubits. 

In the following we define a graph representation of in¬ 
teraction ( |I^ which takes into account this property. Let 
us consider an oriented hypergraph G = {V,E), where 
V = {ij\i < j,j e {1, ■ • ■ ,1V}}. Thus vertices of V rep¬ 
resent pairs of qubits, for instance vertex ij represents 
the pair of qubits i and j. Thus \V\ = (^). If pij^k > 0, 
then the vertices ij,ik, and jk are joined with a F-arc 
whose tail is the vertex ij and whose head is formed by 
vertices ik and jk. We call G the interaction F 2 -graph 
of ( [T7| . An example of the interaction F 2 -graph is de¬ 
picted in Fig. This representation allows us to fully 
describe base F 2 -graphs and the base attractor space by 
the following theorem, which is given without detailed 
proof (see appendix). 

Theorem V.2 Base F 2 -graphs of the considered family 
of interactions are F 2 -graphs with the following property. 
The F 2 G = (V, E) is a base graph if and only if for 
each vertex ij G V and for all k G {1, ■ •. N}\{i,j} exists 
Ik G {I, ■ • ■ N} such that vertices ij and klk are connected 
by a path. The basis of common eigenvectors is given by 
N + 2 vectors \0 n), and Hlj.Af) \i G {1,.. ■,N}} 


with |0,v) o.nd defined by theorem V.l and 


lli.jv) = |0)® ■••(g)|L)(g)---(g)|0). 


24 23 



FIG. 2: An example of an interaction F 2 -graph G correspond- 
ing to the RUO '. This F 2 -graph belongs to the class of 
base graphs of RUO acting on four qubits. This graph 

fulfils the condition given by theorem |V.2| and it thus belongs 
to the set of base graphs. 


All these eigenvectors correspond to the eigenvalue a = 
1 and consequently all base p-attractors solve attractor 
equations © with eigenvalue A = 1. Together with the 
identity operator they form the base attractor space for 
any N > 3. 


Equipped with theorem (V.2), we conclude that in the 
case N > 3 RUO ( |l7| associated with an interaction base 
F 2 -graph generates iterated evolution with the stationary 
asymptotic state Here, 7^0.2) jg projector onto 
the {N + 2)-dimensional space of common eigenvectors. 


C. Reaching consensus within qubit networks with 
cu-interactions 

From a different perspective quantum networks con¬ 
stitute a well designed model of stochastic distributed 
quantum processing. In order to achieve a certain goal, 
individual distant agents of a quantum network inter¬ 
act randomly with their local neighbors according to the 
prescribed interaction graph. One of basic tasks of dis¬ 
tributed processing is reaching consensus among inter¬ 
acting parties. Based on the given situation the desired 
common consensus can take different forms. Here, in the 
context of quantum networks, we follow the definition 
discussed in details in El- Consensus is achieved if the 
total state of the quantum network is permutationally 
invariant, i.e. reordering of qubits of the network does 
not lead to a change of its state. Consequently, it induces 
that all subsystems of the same size are in the same state 
and thus the expectation value of any local observable 
provides the same value irrespective which subsystem is 
being measured. The problem of reaching consensus in 
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quantum networks with swap unitary interaction were 
successfully resolved in [MUTB] . 

In view of Theorems |V.1| and |V.2[ we find out that 
all attractors of quantum networks with a cu-interaction 
base graph are permutationally invariant. Therefore, any 
initial state of such quantum network is asymptotically 
driven into consensus, i.e. quantum network asymptot¬ 
ically reaches consensus. However, we see that the con¬ 
sensus varies for quantum networks with different cu- 
interactions. The resulting structure of asymptotic states 
is determined by a specific information about the initial 
state which is preserved during the evolution. One can 
characterize this memory of initial conditions by the set 
of integrals of motion, i.e. the set of hermitian observ¬ 
ables whose expectation values remain constant during 
the evolution. As the basis of the set of attractors cor¬ 
responding to eigenvalue A = 1 can be always chosen 
hermitian, using Theorems |V.1| and |V.2| one can directly 
construct the whole algebra of integrals of motion of a 
corresponding quantum network. Moreover, if the inter¬ 
action topology of the network is given by a base graph 
the preserved information is uniformly distributed among 
all subsystems and the consensus is established. An in¬ 
teresting question is whether this information can be ex¬ 
tracted by a well designed local measurement solely. The 
following example illustrates that in some cases it is pos¬ 
sible. We assume the observable X = |Ojv)(OAr|, whose 
expectation value Tr[yo(n)A] stays constant during the 
evolution of any quantum network with cu-interaction. 
Further, we consider a quantum network with two-qubit 
cnot-interaction whose qubits interact according to an in¬ 
teraction base graph. We define a local observable acting 
on one qubit 


VI. ASYMPTOTIC EVOLUTION OF 

SIMULTANEOUS TWO AND THREE-QUBIT 
CU-INTERACTION 

Results of the previous section enable us to study 
asymptotics of networks with simultaneous two and 
three-qubit cu-interactions. In such network we assume 
that during each iteration, qubits are allowed to undergo 
either two or three qubit cu-interaction. As this process 
is random, one step of evolution is described by the RUO 
defined as 

$(sim) ^p(2)^(2) _^p(3.1)^(3.1) ^(3.2) ^(3,2) (- 20 ) 

with -I--I-= 1. Using ([^ one can eas¬ 
ily deduce that if all three probabilities and 

p{3,2) nonzero, the corresponding attractor space 
Atr($(®™)) obeys the relation 

Atr($('*™)) = Atr($(^)) n Atr(<l)*^^’^)) n Atr($(^’^^). 

The comparison of base attractor spaces correspond¬ 
ing to three-qubit cu-interactions with the base attractor 
space of two-qubit cu-interactions uncovers several inter¬ 
esting physical implications. Both base attractor spaces 
of three-qubit cu-interactions fulfil relations 

Atr($(^^) C Atr($*^^’*^), i G {1,2}, (21) 

which give us the important attractor space equivalence 



Using Q straightforward algebra reveals equality of ex¬ 
pectation values of both observables 

Tr [p{n)X] = Tr [p,cr], (19) 

which applies to any initial state p(0), any time n S N 
and any reduced one-qubit density operator pi of the 
asymptotic state px . Apparently, mutual two-qubit 
cnot-interactions uniformly distribute information about 
the overlap of the initial state p(0) with the state jO^v) 
among distant constituents of a strongly connected net¬ 
work. This information can be extracted by means of 
local measurements. We emphasize that if the network 
is not strongly connected this property is not guaranteed 
for all initial states p(0) any more. 


d(p(«m)^p(2)) 



FIG. 3: Convergence of the state p^“™^(n) corresponding to 
a network consisted of 4 qubits equipped with simultaneous 
two and three qubit cu-interactions towards the state of the 
same network equipped with purely two-qubit cu-interactions 
We plot the distance d{A,B) — ||A— B\\hs between 
states p^^\n) and for an ensemble of 500 randomly 

chosen initial states, each inserted in operations ( |13[ l and ( |20[ |. 
All relevant graph structures describing the topology of the 
network belong to corresponding sets of base graphs. 












Atr($("™)) = Atr($(2)), 


This allows us to identify the privileged role played by 
two-qubit cu-interactions. Suppose both two and three 
qubit cu-interactions are allowed in the given qubit net¬ 
work. Then, if the interaction graph corresponding to the 
allowed two-qubit cu-interactions is strongly connected, 
the resulting class of possible asymptotic states is not af¬ 
fected by three-qubit interactions at all. After sufficiently 
long time we are not able to recognize, whether observed 
states are results of two-qubit cu-interactions only, or if 
two-qubit cu-interactions are accompanied also by three- 
qubit cu-interactions. Figurej^illustrates this asymptotic 
behavior for 500 randomly chosen initial states. 

Another crucial property is that base attractor spaces 
of considered three-qubit cu-interactions fulfill the rela¬ 
tion 


Atr($(3.i)) n Atr($(3,2)) = Atr($(2)). ( 22 ) 


This relation is displayed schematically in Fig. It 
implies that in the asymptotic regime, three-qubit cu- 
interactions can be indistinguishable from two-qubit cu- 
interactions. The state of a qubit network, in which 
qubits interact exclusively by three-qubit cu-interaction 
can converge towards the asymptotic state of a qubit net¬ 
works with two-qubit cu-interactions. In particular, if 
both interaction F^-graphs belong to corresponding sets 
of base F^-graphs, the possible class of asymptotic states 
is given by the base attractor space of RUO (13). An 


example of this behavior is portrayed for 500 randomly 



FIG. 5: Relation between the base attractor space Atr(<I>l^l) 
of two-qubit cu-interactions (black line) and base attractor 
spaces Atr($*®’*^) of three-qubit interactions with i control 
qubits, i € {1,2} (green and blue plane). 


chosen initial states in Fig. However as simulations 
indicate, this condition is not necessary. Even if one of 
interaction F^-graphs does not belong to corresponding 
sets of base graphs, the asymptotic dynamics can still 
tend towards the one of two-qubit cu-interaction for any 
initial state. 


VII. PROPERTIES OF ASYMPTOTIC STATES 
OF THREE-QUBIT CU-INTERACTIONS 





This section provides additional results concerning 
qubit networks with purely three-qubit cu-interactions. 
The main focus is given to properties of three-qubit cu- 
interactions with one control qubit equipped with a base 
Fi-graph. 


Behavior of qubit networks equipped with interactions 


(14) is according to theorem V.l dependent on the value 
of the parameter (j). Let us consider the parameter (p = ^. 
One can easily see, that in this case interactions (14) pre¬ 
serve the parity of r(z) (defined in theorem V.l). This 
results in the existence of two distinct mixed states In,e 
and In,o belonging to the base attractor space. Further¬ 
more, for ^ = I vectors have a special form. They 
can be written as 


FIG. 4: Convergence of the state p^^\n) corresponding to 
a network consisted of 4 qubits equipped with simultaneous 
three qubit cu-interactions of towards the asymptotic state of 
the same network equipped purely two-qubit cu-interactions 
We plot the distance d{A, B) = \\A — B\\hs between 
states (n) and p^ for an ensemble of 5 00 r andomly chosen 
initial states, each inserted in operation (201 with = 0. 
Both relevant F^-graphs describing the interaction topology 
belong to corresponding sets of base graphs. 


ZGI(N) 

= tJv E I-) ■ 

zei(N) 


It is however more convenient to use a different basis. We 
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thus define 


lx^) = 


V2^-i - 1 


1 ^) 


Z e I{N), 
t{z) — 2fc, 
fe / 0 


= E I")- 


2 G I(N), 
t{z) — 2k — 1 


Together with the vector |Ojv), these form an orthonor¬ 
mal basis of the decoherence-free subspace 
Moreover, this decoherence-free subspace can be di¬ 
vided into two subspaces with the help of projectors 
= |0^) (OatI -h \x%) {x%\ and = \x%) {x%\ 
as 0 Let us denote 

the orthogonal complement to in the subspace of 

Jff with orthonormal basis {|z) |2 e I{N),t{z) = 2k,k G 
MUjO}} and similarly Vq’^^ the orthogonal complement 
to Vq’^^ in the subspace of with orthonormal basis 
{\z) \z G I{N),t{z) = 2k — l,k G M}. The asymptotic 
state can be then written as 




Tr 

V 

P’piO) 

Tr 

' E 



(3.1) 


Tr 

V 


Tr 

' O 



r 


-pCS,!) 
' O 


(23) 


The relation (231 is an analogy of the asymptotic state 
([^. Here, the maximally mixed state living on the or¬ 
thogonal complement of common eigenvectors is splitted 
into two by itself invariant maximally mixed states living 
on subspaces and respectively. 

Another interesting result can be obtained by inspect¬ 
ing the reduced density matrix corresponding to a sin¬ 
gle qubit subsystem. If we define parameters po = 
(OAf|p(0)|Ov) and p± = (v5^|p(0)|(^^), the asymptotic 
state of a single qubit subsystem of a quantum network 
consisting of A^ 3> 1 qubits equipped with three-qubit 
cu-interactions with one control qubit for any </> G (0, tt) 
can be written as 


A3.1).(l) 

r oo 


1 

2 


(Ji + a • (t) , 


with the Bloch vector [3] 


(24) 


S(A :B) 



FIG. 6: The index of correlation of arbitrary two qubits A and 
B in qubit network interacting with two or three qubit cu- 
interactions for different cases. The case po = 0, = p = p_ 

corresponds to the red (blue) line for two-qubit cu-interactions 
(three-qubit cu-interactions with one control qubit). The case 
Po = 1—p, p+ — p = p- corresponds to the green (yellow) line 
for two-qubit interactions (three-qubit cu-interactions with 
one control qubit). 


generally depend on the value of the parameter (j) with 
the exception of the case po = 0. 

Let us now look at entropy related properties of asymp¬ 
totic states two and three qubit cu-interactions. A con¬ 
sequence of the inclusion C with i G {1,2} 

is that for identical initial states p(0) the von Neumann 
entropy S{p) = —Tr[p Inp] [IH3] of the asymptotic states 
Ptc and satisfy the relation S(p^^) > S'(p£’*^). 

Concerning three-qubit cu-interactions with one control 
qubit, there is an additional interesting results concern¬ 
ing the index of correlation 5'(A : B) = S'(A) + S{B) — 
S{A,B) [2] of two arbitrary qubits A and B. The ex¬ 
istence of two orthogonal factorizing eigenstates |(p^) 
results in the possibility of significantly higher correla¬ 
tions between pairs of qubits. This difference between 
two and three qubit cu-interactions is shown in Figfor 
N ^ 1. The maximal value S{A : B) = 1 is achieved for 
three-qubit cu-interactions with one control qubits and 
for value of corresponding parameters p+ = P- = \- 
This is due to the fact that for p+ = p- = ^, the single 
qubit reduced state is according to (24) always the max¬ 
imally mixed state. Also, for any </> G(0, tt) the eigenval¬ 
ues of the two qubit reduced state for p+ = p_ = f are 

{iio,o}. 


VIII. CONCLUSIONS 


a = ((P-r -P-)An(j), 0, po + {p+ -p-)cos(()), 

and (7 = {X,Y,Z) being the vector whose elements are 
Pauli matrices [T]. Eigenvalues of the density matrix 

These eigenvalues thus 


We investigated the asymptotic dynamics of qubit net¬ 
works of arbitrary finite size equipped with two and 
three-qubit cu-interactions for a special class of inter¬ 
action topologies. The main mathematical results are 
presented in the form of two theorems (section which 
characterize asymptotic states of the considered interac- 
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tion topologies. The essential steps of the proof are given 
in appendix. 

For the description of the topology of a given qubit 
network with three-qubit cu-interactions we used the con¬ 
cept of hypergraphs. Although the previously developed 
formalism [TB] of interaction graphs could be adapted for 
considered cases by adding coloring of edges, hypergraphs 
provide a natural way for the description of the topology 
of the interaction. However, in the case of three-qubit 
cu-interaction with two control qubits, there is a signif¬ 
icantly different interpretation of vertices of the interac¬ 
tion graph. This change is necessary for the classification 
of the network topology. The described procedure can be 
extended to general multi-qubit cu-interactions. In this 
respect, the given framework is a practical generalization 
of the interaction studied previously by Novotny et al. 

m- 

Concerning the physics of qubit networks with studied 
interactions, results can be summarized in the following 
way. First, it is shown in section [VT] that enabling three- 
qubit cu-interactions next to two-qubit interactions asso¬ 
ciated with a strongly connected interaction graph leaves 
the possible class of asymptotic states unchanged. This 
is a particularly important result as it implies that quan¬ 
tum networks in which qubits interact primarily via two- 
qubit cu-interactions mask the effect of qubit interactions 
of higher order which thus do not affect the resulting 
asymptotic state. 

Second, relations between base attractor spaces of two 
and three-qubit cu-interactions entail that in the asymp¬ 
totic regime, simultaneous action of both types of three- 
qubit cu-interactions can be indistinguishable from two- 
qubit cu-interactions. This result holds for even more 
general attractor spaces of three-qubit cu-interactions. 

Third, concerning purely three-qubit cu-interactions, 
results given in section |VII| prove that they have dif¬ 
ferent asymptotic properties compared to properties of 
controlled two-qubit interactions. Apart from the above 
mentioned modification of the graph structure used for 
the description of the topology, new attractors emerge in 
both cases. In the case of controlled three-qubit interac¬ 
tion with one control qubit, this leads to the existence 
of two orthogonal factorizing eigenvectors \^%)- Their 
presence in attractor space is responsible for stronger cor¬ 
relations between pairs of qubits than for two-qubit cu- 
interaction. Furthermore, for three-qubit cu-interactions 
with one control qubit, one can distinguish two differ¬ 
ent cases (j> = ^ and (j) ^ ^ with different asymptotic 
properties. The characteristic feature of three-qubit cu- 
interactions with two control qubits is the explicit depen¬ 
dence of the dimension of the base attractor space on the 
size of the quantum network. This property is absent in 
all other considered cases. The additional eigenvectors 
which emerge in this case, have a simple structure 
similar to the eigenvector |0jv). Their existence follows 
directly from the form of the considered interactions. To 
summarize, quantum networks with purely three-qubit 
cu-interactions with a fixed number of control qubits have 


in general a broader class of asymptotic states. 

The studied qubit networks exhibit additional interest¬ 
ing features. If each qubit can ’’see” all other qubits, all 
asymptotic states are permutationally invariant. Hence, 
a qubit network equipped with a base interaction graph 
asymptotically establishes consensus among interacting 
participants. 

Proofs of theorem IV.II and theorem IV.21 indicate the 
existence of general features of asymptotic properties 
of multi-qubit cu-interactions. The number of control 
qubits plays a crucial role in the description of networks 
topology, which can be done analogously to examined 
three-qubit cu-interactions. Increasing the number of 
target qubits can lead to new properties, e.g. the de¬ 
pendence of the form of the base attractor space on the 
value of the parameter (j) or to the existence of additional 
attractors, but it does not change the structure of base 
graphs, which are formed by all strongly connected in¬ 
teraction F-graphs. 

There are several open questions linked to qubit net¬ 
works with random unitary operations. Among them the 
most important question is the connection between the 
probability distribution {pi} appearing in (I) and the 
rate of convergence of the initial state p{0) towards the 
asymptotic state Pao{n). Although several numerical cal¬ 
culations were made, therels no clear mathematical result 
concerning this dependence. Another group of questions 
is linked to consensus formation in qubit networks with 
different interaction topologies. These questions will be 
addressed in future work. 
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Appendix: Some notes about proofs of theorem V.l 
and theorem V.2 

The proof of theorem V.l and theorem V.2 follows con¬ 
cepts developed in [18]. The purpose of this appendix is 
to point out important moments of the proof and main 
differences with proof given in |18j . For a detailed proof 
see [22|. 

As the steps of proofs are similar, we start the proof 
for general three-qubit cu-interaction. The base attractor 
space definitely appears in the case, where every possible 
triplets are present in the set E of an interaction F^-graph 
G = (V, E). Such interaction F^-graph is called the maxi¬ 
mal interaction F^-graph. Elements of computational ba¬ 
sis are coupled with each other by operators C/g. For the 
maximal interaction F-graph we can divide the pairs of 
these elements (which correspond to the elements of den¬ 
sity matrix) to the classes Ai with the property, that pairs 
in each class are all coupled to each other pairs in the 
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FIG. 7: The index graph associated with the interaction 
F-graph and the corresponding index graph 



In the case of three-qubit cu-interactions with one con¬ 
trol qubit a single operator Ue couples up to four elements 
of the computational basis depending on the particular 
element and the value of the parameter cj). For a given 
index graph corresponding to the three-qubit cu- 

interaction with one control qubit one can find a cer¬ 
tain index graph corresponding to the two-qubit cu- 
interaction which is a subgraph of the former index graph 
(apart from the case cj) = ^). An example of this con¬ 
struction is displayed in Fig. [Tj 

Afterwards, we prove that one of these index graphs 
is two-connected if and only if the other one is two- 
connected thus completing the proof for the case (p ^ 
|. For the case ^ | the index graph always 

has at least two components. We find that the in¬ 
dex graph then corresponds to the base Fi-graph if 
and only if it has exactly two components which are 
two-connected. We can prove the theorem for this 
case by taking the generalization of the so-called star 
graph (which we call star Fi-graph) [^. The star 
Fi-graph is the Fi-graph G = {V,E) with V = 
andii; = {(l,2j)|j € {3,... N}}[j{{2,lj)\j G 


same class and they are not coupled to any pairs in other 
classes. To each of these classes corresponds a single lin¬ 
early independent attractor which is then found directly 
by solving attractor equations. Couplings between these 
pairs can be visualized by the so-called index graph [TH] . 
The index graph is an undirected colored graph, whose 
vertices represent elements of computational basis. In 
its definition we omit all elements of computational basis 
which correspond to the classes Aj which are formed by 
a single element. Two vertices i and j of a given index 
graph are connected by an edge of color e, if the elements 
|z) and |j) of the computational basis are coupled by the 
operator Ue- 

Now depending on the particular type of interaction 
we need certain subgraph of the index graph to be two- 
connected to obtain the base F^-graph. Next, analo¬ 
gously to the procedure used in [18] we prove by induction 
that we can eliminate hyperedges from the interaction F^- 
graph as long as the newly emerged interaction F^-graph 
fulfils the statement of the corresponding theorem. 



FIG. 9: The set Q 4 , corresponding to the star F 2 -graph on 4 
vertices. The star F 2 -graph is used in the proof of theorem 

ED 
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{3,... N}} lJ{(j, 12)|j e {3,... N}}. An example of star 
Fi-graph on 4 vertices is given in Fig. Equipped with 
this Fi-graph, we follow the same induction steps as in 
the proof of the form of base graphs of two-qubit inter¬ 
actions HD- 

The case of three-qubit cu-interactions with two con¬ 
trol qubits is treated similarly. However, the generaliza¬ 
tion of the star graph (which is called star F 2 -graph) 
needs to be constructed in a different way. First, for an 
interaction F 2 -graph G = {V, E) we define a set of ordi¬ 
nary directed graphs Q = {Gi = {Vi,Ei)}. Here Gi is 


an ordinary directed graph with the set of vertices Vi = 
{j\j G {1, ■ • ■ ^}\{0} and edges {j, k) e Ei k) G 

E. We define star F 2 -graph as follows: G = {V,E) is a 
star F 2 -graph, if the corresponding set of graphs Q con¬ 
tains two graphs, which are ordinary star graphs and the 
rest of graphs contained in Q are isomorphic with each 
other. Figure shows the set corresponding to star 
F 2 -graph on 4 vertices. After constructing the general¬ 
ized star graph the proof follows the same steps as are 
used in [18]. 


[1] M. A. Nielsen and I. L. Chuang. Quantum Computation 
and Quantum Informaion. Cambridge University Press, 
2010 . 

[2] S. M. Barnett. Quantum Information. Oxford University 
Press, 2009. 

[3] B. Schumacher and M. Westmoreland. Quantum Pro¬ 
cesses Systems, & Information. Cambridge University 
Press, 2010. 

[4] M. Ballerini et al, Proc. Nat. Acad. Sci. USA 105, 1232- 
1237, 2007. 

[5] C. Huepe et al, New J. Phys. 13, 073022, 2011. 

[6] H. J. Kimble, Nature 453, 1023-1030, 2008. 

[7] J. Gough and M. R. James, IEEE Trans. Autom. Control 
54, 2530-2544, 2009. 

[8] X. Wang, P. Pemberton-Ross, and S. G. Schirmer, IEEE 
Trans. Autom. Control 57, 1945-1956, 2012. 

[9] I. Cirac, P. Zoller, H.J. Kimble, and H. Mabuchi, Phys. 
Rev. Lett. 78, 3221, 1997. 

llOl A. Acln, I. Cirac, and M. Lewenstein, Nature Physics 3, 
256-259, 2007. 

[11] S. Ritter et al., Nature 484, 195-200, 2012. 


[12] J. Novotny, G. Alber, and I, Jex, Phys. Rev. Lett 107, 
090501, 2011. 

[13] F. Jin et al, Phys. Rev. A 87, 022117, 2013. 


L. Mazzarella, 
arXiv: 1303.4077, 

A. 

Sarlette, 

and 

F. 

Ticozzi, 

L. Mazzarella, 

A. 

Sarlette, 

and 

F. 

Ticozzi, 


arXiv:1403.3582vl. 

[16] G. Shi, D. Dong, and K. H. Johansson, 
arXiv:1403.6387v3. 

[17] J. Novotny, G. Alber, and 1. Jex, Central European J. 
Phys 8, 1001, 2010. 

[18] J. Novotny, G. Alber, and 1. Jex, New J. Phys. 13, 
053052, 2011. 

[19] G. Gallo, G. Longo, S. Nguyen, and S. Pallottino, Dis¬ 
crete App. Math. 42, 177-201, 1993. 

]20] B. Kollar, J. Novotny, T. Kiss, and 1. Jex, Eur. Phys. J. 
Plus 129, 103, 2014. 

[21] R. Diestel, Graph Theory, Springer, 2000. 

[22] J. Maryska, Ing., Quantum networks with three particle 
interactions, Master’s thesis, FNSPE CTU Prague, 2013. 



